Gapless layered three-dimensional fractional quantum Hall states 
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Using the parton construction, we build a three-dimensional (3D) multilayer fractional quantum 
Hall state with average filling v = 1/3 per layer that is qualitatively distinct from a stacking of 
weakly coupled Laughlin states. The state supports gapped charge e/3 fermionic quasiparticles that 
can propagate both within and between the layers, in contrast to the quasiparticles in a multilayer 
Laughlin state which are confined within each layer. Moreover, the state has gapless neutral collec- 
tive modes, a manifestation of an emergent "photon" , which is minimally coupled to the fermionic 
quasiparticles. The surface sheath of the multilayer state resembles a chiral analog of the Halperin- 
Lee-Read state, which is protected against gap forming instabilities by the topological character 
of the bulk 3D phase. We propose that this state might be present in multilayer systems in the 
"intermediate tunneling regime", where the interlayer tunneling strength is on the same order as 
the Coulomb energy scale. We also find that the parton construction leads to a candidate state for 
a bilayer v = 1/3 system in the intermediate tunneling regime. The candidate state is distinct from 
both a bilayer of v — 1/3 Laughlin states and the single layer v — 2/3 state, but is nonetheless a 
fully gapped fractional quantum Hall state with charge e/3 anyonic quasiparticles. 



I. INTRODUCTION 

While the fractional quantum Hall effect is firmly 
rooted in two dimensions (2D), anisotropic 3D electron 
systems - such as multilayer systems in a perpendicu- 
lar magnetic field - can exhibit fractional quantum Hall 
(FQH) states, at least in principle. The simplest example 
of such a 3D multilayer state is a stacking of N decoupled 
Laughlin states. [H,0] At the next level of complexity, one 
can construct states with interlayer correlations, such as 
the (3, 3, 1) bilayer state. [H More generally, multicom- 
ponent Chern-Simons theory allows one to construct a 
myriad of TV-layer analogues of the (3, 3, 1) state. [4| 

These states are quite general but they suffer from a 
limitation: they all have a fixed number of electrons in 
each layer. This restriction could be problematic for de- 
scribing certain multilayer systems, especially those with 
appreciable interlayer tunneling. Therefore, alternative 
constructions of 3D multilayer FQH states are desirable 
theoretically. 

On the experimental side, a number of experiments on 
3D semiconductor multilayers have explored the behavior 
of stacked integer quantum Hall states [E B 0] > includ- 
ing the novel vertical transport due to the conducting 
surface sheath. 0, i, d [TO, [u, E 0, 03 More recently, 
experiments on bismuth crystals in high magnetic fields 
have revealed intriguing anomalies in the ultra quantum 
limit - the limit where the magnetic field is sufficiently 
large that only the lowest Landau level is (partially) oc- 
cupied. [l5[ It has been suggested that a novel 3D frac- 
tional quantum Hall type state might be present. While 
these are not layered materials, strong electron correla- 
tions could drive a transition wherein the electrons spon- 
taneously form a weak layered structure, as suggested in 
recent work. [1 , [ItJ Bulk graphite also reveals transport 
anomalies in the ultra quantum limit [la , 1 19L |20| , which 
have been attributed to a charge density wave transition 
in this layered material. The observed quantum Hall ef- 
fect in graphene [2l|, [22|, [H, [24[ , and future prospects for 



graphene multilayers, provides further impetus to explore 
3D layered quantum Hall phenomena. 

Motivated by these experiments, and the lack of pre- 
vious theoretical exploration, we revisit the behavior of 
multilayer systems in the fractional quantum Hall regime. 
Generally, we are interested in addressing the following 
class of questions: what is the fate of a weakly coupled 
stacking of 2D fractional quantum Hall states when the 
interlayer electron tunneling becomes strong enough to 
close the quantum Hall gap? In particular, are new 
fractional quantum Hall type states possible when the 
Coulomb interaction is comparable to the interlayer tun- 
neling strength? We believe that this intermediate tun- 
neling strength regime is both experimentally accessible 
and theoretically novel. 

Answering these questions definitively for a specific mi- 
croscopic model is quite challenging and likely requires 
extensive numerical calculation. Here, we are less ambi- 
tious. Our goal is simply to find candidate ground states 
for the intermediate tunneling regime. This is already a 
non-trivial problem since, as we mentioned earlier, most 
multilayer states - such as those obtained from Chern- 
Simons mean field theory - have a fixed number of elec- 
trons in each layer and hence are unnatural unless the 
interlayer tunneling is weak. 

In this paper, we construct a candidate ground state 
for the simplest possible multilayer system: spinless (or 
spin polarized) electrons with an average filling of v = 
1/3 per layer. We speculate that the candidate state may 
be realized at intermediate tunneling strength. However, 
our arguments for the candidate state are indirect, as we 
do not make any detailed analysis of energetics. 

We construct our candidate state using a slave-particle 
gauge theory approach known as the "parton construc- 
tion." [l^,[2y| The basic idea of the parton construction 
is to write the electron creation operator as a product of 
several (in our case, 3) fermionic parton creation opera- 
tors. By choosing different mean field parton states, and 
including (gauge) fluctuations, one can construct differ- 
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ent FQH states. The advantage of the parton construc- 
tion is that it naturally leads to states with electron num- 
ber fluctuations in each layer. Thus, parton FQH states 
may be particularly natural in the intermediate tunnel- 
ing regime. In addition, the particular state we construct 
has the interesting property that one can tune from it to 
a decoupled multilayer Laughlin state by changing a sin- 
gle coupling constant in the parton gauge theory. Given 
that the multilayer Laughlin state is likely realized at 
weak interlayer tunneling, this is an additional reason 
we consider our state to be a candidate ground state for 
intermediate interlayer tunneling. 

We analyze the candidate state for both a hnite num- 
ber of layers N, and for the 3D limit, N — > oo. For a 
finite number of layers N, we find that the state is a fully 
gapped FQH state (when N = 1, it is simply the Laugh- 
lin state). The quasiparticle excitations are anyonic and 
carry charge e/3. We find that the quasiparticle exci- 
tations are described by a if-matrix [271 ] of dimension 
(3N — 2) x (3JV — 2), along with a charge vector of length 
(37V -2). 

In the 3D limit, the candidate state exhibits more 
unusual physics. It supports two types of excitations: 
gapped charge e/3 fermionic quasiparticle excitations 
and gapless, electrically neutral collective modes. The 
e/3 fermionic quasiparticle excitations (which are essen- 
tially the "deconfined" partons) are truly 3D quasipar- 
ticles and can move freely between layers. This should 
be contrasted with the e/3 excitations in the multilayer 
Laughlin state which are confined to individual layers. 
(In this sense, our candidate state is "more 3D" then 
the multilayer Laughlin state). As for the gapless elec- 
trically neutral collective modes, these are the emergent 
U(l) gauge bosons or "photons" which originate from 
fluctuations about the mean field parton state. Unlike 
Maxwell photons, these excitations have only one polar- 
ization state and have an anisotropic dispersion of the 
form W 2 . ~ + k z (for layers oriented in the xy plane) . 
The e/3 fermionic quasiparticles are minimally coupled 
to these "photon" modes and thus have long range inter- 
actions. 

The edge physics of the TV-layer and 3D systems is also 
interesting. In the case of a finite number of layers N, the 
edge theory is a chiral boson conformal field theory with 
37V — 2 chiral modes. The edge Lagrangian can be read 
off from the bulk if-matrix using the standard formalism. 
[13] The edge (or surface) physics in the 3D limit is more 
complex. For a 3D system with layers in the xy plane and 
a boundary in the xz plane, the three flavors of dispers- 
ing edge modes form a "sheath" of chiral fermions. At 
the mean field level these fermions are non-interacting, 
but with fluctuations included are minimally coupled to 
the gapless bulk "photons" . The surface sheath, then, re- 
sembles a chiral analog of the Halperin-Lee-Read state, 
[28| which is protected against gap-forming instabilities 
by the topological character of the bulk phase. 

This paper is organized as follows. In section |TT] 
we speculate about phase diagrams of multilayer v = 
1/3 systems and we describe the intermediate tunneling 
regime in more detail. In section IIIII we construct our 
candidate state. In sections HVl - I VIA we analyze the bulk 



physics of the candidate state for single layer, bilayer, TV- 
layer, and 3D systems. In section [Villi we investigate the 
crossover between 2D and 3D physics in systems with a 
large but finite number of layers. In section IIXI we de- 
scribe the relationship between the candidate state and 
the multilayer Laughlin state. Finally, in section [X] we 
analyze the edge physics of the candidate state for N- 
layer and 3D systems. 

II. MODEL AND POSSIBLE PHASE 
DIAGRAMS 

In this section, we discuss the physics of iV-layer v = 
1/3 FQH systems in more detail. We speculate about 
possible phase diagrams and we explain what the inter- 
mediate tunneling regime is, and why it is interesting. 

Consider a geometry where the layers are oriented in 
the xy plane, and neighboring layers are spaced a dis- 
tance a in the z direction. There are four energy scales 
in the problem - two intralayer and two interlayer scales. 
The intralayer energy scales are the cyclotron energy 
huj c = heB/m, and the characteristic intralayer Coulomb 
energy scale E c = e 2 /Zg. The interlayer scales are the 
interlayer Coulomb energy e 2 /a and the interlayer tun- 
neling strength t z . In the following discussion, we will 
focus on the regime where (1) hio c is much larger than 
any of the other energy scales, and (2) a is comparable 
to, but larger than lg so that the interlayer Coulomb 
energy scale is smaller than, but on the same order as 
the intralayer energy E c . In this regime, there is only 
one dimensionless parameter in the problem: the ratio 
9 = t z /E c . 

Let us think about the phase diagram as we vary the 
dimensionless ratio g = t z /E c . To begin, suppose N = 2. 
When g = 0, there is no interlayer tunneling and we 
expect that the ground state is given by two decoupled 
Laughlin v — 1/3 states, with perhaps small quantitative 
modifications due to the interlayer Coulomb interaction. 
Since the Laughlin state is gapped, it will be stable to 
small interlayer hopping, i.e. g <C 1. In the opposite 
limit, with very large interlayer tunneling j » 1, all of 
the electrons will be in the "bonding band" - the band 
consisting of symmetric combinations of Landau orbitals 
in the two layers. The system is thus an effective single 
layer system at filling v = 2/3. The weak Coulomb inter- 
actions will presumably lead to an abelian iy — 2/3 state 
with gap of order E c . 

Now, consider the regime g ~ 1. Starting from the 
decoupled limit, when g is increased one expects that the 
quasiparticle gap in each layer will diminish and presum- 
ably be driven to zero at some critical value g\. On the 
other hand, when g is brought down from large values, 
the gap of the "single layer" v = 2/3 state will dimin- 
ish (due to Landau level mixing into the "antibonding 
band") and be driven to zero at some value, 172- Gener- 
ally, there is no reason to expect that g\ — g2, although 
it is possible that there is a direct first order transition 
between the 1/3 + 1/3 decoupled phase and the 2/3 "sin- 
gle" layer state. If gn > gi, there will be a third FQH 
phase (or phases) for g\ < g < g^. This potential phase 
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is the "intermediate coupling" phase we consider in this 
paper. 

One can imagine a similar scenario for the 3D limit, 
N — > oo. Again, at weak interlayer tunneling, j « 1, 
the stack of Laughlin states is stable and can be readily 
analyzed. [l[ On the other hand, the nature of the strong 
tunneling phase with g 3> 1 is non-trivial. In the non- 
interacting limit, the lowest Landau level in each layer 
will form a band that disperses in the z— direction, with 
energy —t z cos(k z a). The noninteracting ground state 
consists of completely filling all of the lowest Landau 
level band states with \k z \ < kp = ■n/3a. This describes 
a gapless state. In a Wannier basis of orthonormalized 
lowest Landau level wavefunctions with guiding centers 
sitting on the sites of a regular 2D lattice (say triangular), 
one can view the system as an array of one-dimensional 
non-interacting electron systems. Gapless particle-hole 
excitations exist across the two Fermi points in each of 
the ID "wires". This noninteracting state is likely to 
be unstable in the presence of arbitrarily weak Coulomb 
interactions, E c ^ 0, due to the nested 2kp backscat- 
tering interactions between nearby "wires" . The sim- 
plest scenario would be the development of a fully gapped 
Q = 2kp = 2ir/(3a) charge density wave (CDW) state, 
which corresponds to a tripling of the unit cell along the 
z— axis. Naively one would expect the CDW gap to scale 
as A C dw ~ t z exp(-const. ■ g). [H,[30| This CDW state 
can be loosely thought of as an effective system with 7V/3 
layers, each one at filling v e ff = 1- 

As in the bilayer case, increasing g from small val- 
ues will presumably close the Laughlin quasiparticle gap 
in each of the layers, destroying the decoupled phase at 
some <7i. Similarly, when g is decreased from very large 
values down to order one the CDW state will become 
disfavored due to the increasing intralayer Coulomb re- 
pulsion. One expects the CDW state to be destroyed for 
some g < g2- As for N — 2, it is possible that there is a 
third phase for g\ < g < g 2 - an "intermediate tunneling 
phase." 

The possibility of such an intermediate tunneling phase 
for either the finite N case or the 3D limit N — > 00, 
is the starting point for this paper. One reason we 
feel it is a particularly interesting possibility is that the 
interlayer tunneling and the Coulomb interaction are 
both of paramount importance in such a putative phase. 
This poses a theoretical challenge since the obvious FQH 
states - such as those constructed from Chern-Simons 
mean field theory - have a fixed number of electrons in 
each layer and are therefore unnatural except for very 
weak interlayer tunneling. On the other hand, if one tries 
to construct a state by treating the A^-layer system as an 
effective single layer system at filling N/3, the result is 
unnatural except for very strong interlayer tunneling. In 
this paper, we use a different approach - a slave particle 
construction with fermionic "partons" - to build a candi- 
date state that overcomes these difficulties. 



III. PARTON CONSTRUCTION 

Our candidate state can be described most naturally 
using the parton construction. [2^, [2(| Let us first de- 
scribe the construction in the case of the single layer sys- 
tem; we will then generalize to multiple layers. 

Our starting point is the single layer electron Hamil- 
tonian. As it will be convenient in what follows, we reg- 
ularize this Hamiltonian, replacing the 2D continuum by 
a square lattice. We take the flux through each plaque- 
tte in the lattice to be 2ir/M, the electron density to be 
1/3M, and we consider the limit M — > 00. In this limit, 
the lattice model behaves like a 2D continuum with elec- 
trons at filling fraction v = 1/3. 

The regularized electron Hamiltonian can be written 

as 

H = — \J(fcJ.e lA!C,< Cx+xi + h.c.) + interactions (1) 

xi 

where A x .i is a lattice gauge field with AiA^ - &-2A x .i — 
2n/M, and A is periodic with unit cell of size M. Here, 
Aj denotes a lattice derivative in the Xi direction: Aif x = 

fx+Xi fx- 

In the parton construction, we think of the electron as 
a composite of 3 fermionic partons dP , p = 1,2,3: 

c = d 1 d 2 d 3 (2) 

We then substitute the expression for c into this Hamil- 
tonian and expand around a saddle point. The result is a 
non-interacting mean-field Hamiltonian for the partons. 
Many different saddle points can be stabilized depending 
on the details of the interactions. Different saddle points 
correspond to different FQH states. 

Here, we consider a particular saddle point. The saddle 
point we are interested in is associated with the mean 
field Hamiltonian 

Hmf = -J2(t P (d^e lA -d p x+it +h.c.) (3) 

xip 

where A x ,i is a lattice gauge field with AiA E ,2~A 2 A r ,i = 
27r/3M, and A X} i is periodic with unit cell of size 3M. 
Notice that the flux 2tt/3M is exactly the right size so 
that the partons are at filling v = 1 (the partons, like the 
electrons, are at density 1/3M). We have also assumed 
that the hopping amplitudes t p are different for the three 
species of partons. 

What is the physics of this parton state? At the 
mean field level, the parton state is a gapped state with 
fermionic excitations. However, this mean field result is 
not quite correct, since we have not taken into account 
the effect of fluctuations about the saddle point. These 
fluctuations are described by fluctuations in the hopping 
amplitudes t p of the form t p — > t p e l6p with 61+62+03 = 0. 
We can parameterize them in terms of two U(l) gauge 
fields A q , q — 1, 2 by setting 9 P = Q pq A q where q = 1,2, 
and Q p i = (1,0, -1), Q p2 = (0, 1, -1). The effect of fluc- 
tuations is thus to couple the partons to two U(l) gauge 
fields A q . (Note that the structure of the gauge fluctua- 
tions is closely tied to the symmetries of the saddle point. 
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For example, at the symmetric saddle point ti = t% = t% 
the fluctuations are described by SU(3) gauge fluctua- 
tions t p — > tp'IJP rather than the U(l) x U(l) fluctua- 
tions present here [3l|). Including these fluctuations, our 
Hamiltonian is given by 

H = H t +H A (4) 

where H t describes the parton hopping and H a describes 
the gauge field dynamics: 

H t = -Y^(t P (<) t e i ^ A h +iA ^dl +ii +h.c.) (5) 

xip 

Ha = ^f(^,J 2 -E Jcos ( A i<2-A2< 1 )(6) 

xiq xq 

The parton state we are interested in is described by the 
above Hamiltonian in the weak gauge fluctuation regime 
- e.g. g < J, t p . 

Generalizing this construction to the N layer case is 
straightforward. In this case, the mean field parton 
Hamiltonian is given by 

xzip 

- Y / M dP ^<( Z+ D + h - c -) (7) 

xzp 

where z is the layer index and t p j_ , t p3 are the intralayer 
and interlayer hopping amplitudes. 

Including the U(l) x U (1) gauge fluctuations, we arrive 

at 

JV N—l 

H = J2( H ** + H * A ) + E ( H '('+i)t + h *(*+i)a) ( 8 ) 

z=l z=l 

where H zt , H z a describe the intralayer hopping and 
gauge field terms, 

xip 

(9) 

h *a = E y (*w-) 2 - E j -l cos ( A ^L, 2 - a 2 ^L,i) 

xig xq 

(10) 

and -ff z ( z +i)t, -ffzrz+i)^ describe the interlayer hopping 
and gauge field terms, 

xp 

h z ( z+ i)a = E y ( E xz,3) 2 

xq 

£ J 3 cos(A.ML, 3 " ^ + ( n ) 

Again, we assume that the gauge fluctuations are weak - 
.93 < J3,t p 3,tpX and g x < J±,t p3 ,t p± . In the following 
sections, we analyze the physics of this state. 



IV. SINGLE LAYER 

We begin with the simplest case: the single layer par- 
ton state. We rederive the well-known result that the 
sing le layer parton state is precisely the Laughlin state. 

To understand the properties of the single layer parton 
state, we need to analyze the low energy physics of the 
Hamiltonian One way to do this is to introduce U{1) 
gauge fields to describe the parton number currents: 

f X = ^e^d^al (12) 

The low energy effective theory for the parton hopping 
terms H t can then be written as 

L = T~ E ^"dyP^l + minimal coupling to A q (13) 

p 

Including the minimal coupling to A q gives 

L = ^ E + ^ v Q vq A\d,al (14) 

v 

Adding the gauge field terms Ha, expanding the cosines 
to quadratic order, and taking the continuum limit, we 
arrive at the low energy effective theory 

L = ^E eV XM + ^ AA "U^M 

+ E^(^-^o) 2 

iq " 

i 

where / is the lattice spacing. The last two terms are 
irrelevant to the low energy physics since integrating out 
the field produces a Chern-Simons term for A q (which 
has one less derivative then the above Maxwell terms). 
Dropping these terms and integrating out A leaves us 
with 

L =^J2^ v < a ^< ( 16 ) 

p 

together with the constraints c^aj, = d^a* — d^a^. Let- 
ting a v — a* = = we get 

L = ^-e^ u 3a x d^ (17) 

47T 

If we include the coupling to the physical electromagnetic 
gauge field Aem - assigning electric charges ei,e2,e3 to 
the partons with e\ + e-2 + e 3 = e - we find 

L = J_ e AM- 3aAS av + JL e ^ AEM xd aij (18 ) 

(irrespective of the values of ei,e 2 ,e3). This is the low 
energy effective theory for the Laughlin state. We con- 
clude that the single layer parton state is in the same 
universality class (e.g. quantum phase) as the Laughlin 
state. 
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V. BILAYER 

In this section we analyze the parton state in the next 
simplest case: a bilayer. In this case, the parton Hamil- 
tonian (JSJ) reduces to 

2 

H = J2(Hzt + H zA ) + (H 12t + H 12A ) (19) 

2=1 

where H zt , H z a describe the intralayer hopping and 
gauge field terms, 

H zt = - Efe-L^,)^" 3 '^"^-^, + h.c.) 

xip 

(20) 

H zA = Y&kif - E J -L cos ( A ^L,2 " A^L.r) 

(21) 

and Hi2t, H\ia describe the interlayer hopping and gauge 
field terms, 

Hnt = -Y,M€i) ] j Q * qAll ' 3 <2 + h.c.) (22) 

xp 

H 12A = Efffi,s) 2 

xq 

- ^JscosCAM^-^m + ^v) (23) 



xtq 



As before, we can derive the low energy physics of this 
Hamiltonian by introducing U(l) gauge fields a|„ to de- 
scribe the parton number currents in each layer: 



Jz 2^ u H a zv 



(24) 



Expanding the cosines in the gauge field terms to 
quadratic order, and putting everything together, we ar- 
rive at the effective theory 



zp 



+ ( do Al 3 -Al + Al f 

Zg s a 

- Y^ J i{d i Al,-A\. + Al.f 



(25) 



where a is the layer spacing. (As in the single layer case, 
we have dropped the intralayer Maxwell terms as they 
are irrelevant to the low energy physics). To proceed 
further, we choose the gauge A\ 3 = 0, and define new 
fields A q± = A\ ± A\. Expressing the Lagrangian in 
terms of these fields gives 

J — — \^ fVn 1 " ri nP 

4^ 2^ a z,\°H a z,v 

zp 



-\2 



(26) 



As in the single layer case, the final step is to integrate 
out the gauge fields A q± . Integrating out A q+ generates 
the constraints J2 Z <^ a z> = J2 Z dn a l,u = J2 Z ^ a z,u'i in- 
tegrating out A q ~ generates a Maxwell term for a p which 
is irrelevant to the low energy physics due to the presence 
of the Chern-Simons term. 

We thus arrive at the Lagrangian 



zp 



(27) 



l z,v 



together with the constraints J2 Z ^^ a l v = J2 Z ®v a z 

^M a z,y There are four independent gauge fields left 

which we can parameterize by a 1 = a\,a 2 — a 2 , a 3 = 

a?, a 4 = Ez°z = 
variables, we have 



Sz a z = J2z a z- I n terms of these 



L = 



^\pv 



in 



Kua[d, 



(28) 



where 



K 




(29) 



Including the coupling to the physical electromagnetic 
gauge field Aem, assigning charges e\,e2,e$ to the par- 
tons with ei + e2 + e3 = e, we find 

L = Le^RuaidX + ^-e x ^t I A EM . x d ti ai (30) 

47T Z7T 

where t T — (0,0,0,1) (irrespective of the values of 
ei,e 2 ,e 3 ). 

The parton state is completely specified by the above 
.ftf-matrix and charge vector t (or more accurately, the 
universal properties of this quantum state are completely 
specified). We now analyze the basic properties of this 
state. 

We begin with the quasiparticle statistics. According 
to the if-matrix formalism, the quasiparticle excitations 
can be labeled by integer vectors I. The exchange statis- 
tics of a quasiparticle I is given byd ex =tt(I t K-H). The 
mutual statistics of two quasiparticles 1,1' - e.g. the phase 
associated with braiding one particle around another - is 
given by 9 mut = 2n(l T R-H'). 

In principle, these formulas completely specify the 
quasiparticle statistics of the parton state. However, 
it is convenient to describe the statistics of the parton 
state in a more concise way. A general way to do this 
is to find a subset of quasiparticles with the property 
that one can generate all topologically distinct quasi- 
particles by taking composites of these basic quasipar- 
ticles. One can then describe the complete quasiparticle 
statistics by specifying the statistics of this generating 
subset of quasiparticles. For the above state, the three 
parton excitations corresponding to h = (1,0, 0,0), li — 
(0, 1,0,0),Z 3 = (0,0,1,0) generate all the others. (One 
way to see this is to note that the excitation (0, 0, 0, 1) 
is topologically identical to (2,0,0,0)). Simple algebra 
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shows that the three parton excitations have exchange 
statistics 



9-7T 



and mutual statistics 



%) = 



(1 + 35 pp >)iT 



(31) 



(32) 



This gives a complete description of the quasiparticle 
statistics of the bilayer parton state. The electric charges 
of the quasiparticles are also easy to obtain. Again, it 
suffices to specify the parton charges, which are given by 



q p = e ■ (t T K-H p ) = - 



(33) 



Now that we have computed these properties, we can 
see that the bilayer parton state is distinct from a bi- 
layer of decoupled v = 1/3 Laughlin states as well as the 
conventional v — 2/3 state. Indeed, one can distinguish 
the states by noting that the e/3 excitation in the parton 
state has a statistical angle 27r/3, while the e/3 excitation 
in the other two states has an angle 7r/3. 

One can also distinguish the states by their ground 
state degeneracy on a torus. This quantity is particularly 
easy to measure in numerical calculations. The ground 
state degeneracy for the bilayer parton state is just the 
determinant of K which is 12. On the other hand, the 
degeneracy of a bilayer of decoupled v = 1/3 Laughlin 
states is 9, and the degeneracy of the v = 2/3 state is 3. 

A final way to distinguish the states is via their ther- 
mal Hall conductances. Recall that gaph chiral boson 
edge mode gives a contribution of ±^f-T to the ther- 
mal Hall conductance, with the sign determined by the 
chirality of the mode. Thus, the thermal Hall conduc- 
tance can be computed by counting the number of pos- 
itive and negative eigenvalues of K. In the case of the 
bilayer parton state, there are four positive eigenvalues 

so the thermal Hall conductance is 4 (in units of 3h B T). 
On the other hand, the thermal Hall conductance for the 
bilayer of Laughlin states is 2 and the thermal Hall con- 
ductance for the v = 2/3 state is 0. 



VI. iV-LAYER SYSTEM 

The bilayer results can be easily generalized to the A- 
layer case. For general N, one finds a A-matrix of di- 
mension (3 A — 2) x (37V — 2). The result is shown below 
for the case A = 3: 



K = 



i 2 
1 











1\ 



-1 
-1 
-1 
-1 
-1 



(34) 



\-l -1 -1 -1 -1-1 3 / 



The corresponding charge vector is t T = 
(0,0,0,0,0,0,1). The generalization to arbitrary 



A is clear: along the diagonals there are three 
(A - 1) x (N - 1) blocks of the form 1 + while the 
last row and column is made up of — l's with a 3 in the 
bottom right-hand corner. 

As before, the A-matrix and charge vector deter- 
mine all the universal properties of the FQH state, such 
as the quasiparticle statistics and charges. Also, just 
as before, one can summarize the quasiparticle charges 
and statistics more concisely by specifying the statistics 
and charges of the three parton species (which corre- 
spond to h = (1,0,0,0,0,0,0), h = (0,0,1,0,0,0,0), 
h = (0, 0, 0, 0, 1, 0, 0) in the A = 3 case). 

One finds that the parton excitations have exchange 
statistics 



mutual statistics 



and charge 



_ (3A - 2)tt 
p ~ 3A 



(2 + (6A-6)W 
3A 



* = 3 



VII. 3D LIMIT 



(35) 



(36) 



(37) 



In this section, we analyze the parton construction in 
the 3-D limit, A — > oo. Recall that the parton Hamilto- 
nian is given by 

H = *^2(H zt + H zA + ff z ( z+ i)( + H z ( z+ i) A ) (38) 



where 



xip 

(39) 



and 



xiq 



xq 



and 



(40) 



H 



z(z + l)t 



Efe3(rfL) t ^ AL ' 3 < (z+1) + ^ C -) 

xp 
93, 

2 



H Z(Z+1)A = Ef(C, 3 ) 2 (41) 

xq 

- E J 3Cos(A^ z , 3 -AL, + ^ ( , +1) ,) 



xiq 



As usual, we derive a low energy effective theory by in- 
troducing U(l) gauge fields to describe the parton 
number currents in each layer: 



J z ~ 2ir 



(42) 
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Expanding the cosines in the gauge field terms to 
quadratic order, and putting everything together, we ar- 
rive at the effective Lagrangian 



L = L n 



L A 



where 



L a =— y 

zp 



(43) 



(44) 



and 



La = 



ziq 



zq 



z<j y 

E Tr(^^z,3 ~ A l.i + -^z+i.j) 2 



and 



= -E< 



Q 



pg^ ± z,x u ^ u z : i> 



(45) 



(46) 



To proceed further, we integrate out the gauge fields dP z 
corresponding to the parton excitations. Note that this 
is different from the approach we took in the single layer 
and bilayer cases where we integrated out the gauge fields 
A q instead. We could have used this approach in those 
cases as well. The advantage of this approach is that 
it leads to a simpler description of the bulk low energy 
physics: the low energy effective theory for the N layer 
case is simply a 2 x 2 Chern-Simons theory coupled to 
fcrmionic partons (instead of a (37V — 2) x (3N — 2) Chern- 
Simons theory coupled to bosons). The disadvantage is 
that the edge physics cannot be easily read off from the 
bulk effective theory. Here our primary interest is in the 
bulk physics - thus we choose to integrate out the fields 

Integrating out the a field (e.g. the partons) produces 
a Chern-Simons term for the A field. We can then drop 
the J± and g± Maxwell terms, as they are irrelevant at 
long distances. The resulting Lagrangian is given by 

LA,eff = ^J2^K qq ,Al x d,At 

Z 

+ E^( 9 °^,3-^,0 + ^ + l,0) 2 



The full low energy effective theory is described by 
fcrmionic partons minimally coupled to this gauge theory 
with gauge charges Q pq : 



L> Impart 

+ L A,eff 



(49) 



where 



-'part 



J2(d p z )Hido + Q Pq A%)d» z 

pz 

E^t^^-^^m) 2 ^ 

ipz " 

Y J h^e^ A U<P z+1 +h.c. (50) 



(Here, m p ± — l/(2t p ±l 2 )). We now analyze the physics 
of this low energy effective theory. We begin with the 
excitations in the bulk. There are two types of exci- 
tations: gapped parton excitations described by dP and 
gapless gauge boson excitations described by the above 
U{1) x U(l) gauge theory. 

Let us try to understand the gapless gauge boson exci- 
tations in greater detail. We can derive the dispersion re- 
lation for these gapless modes by going to Fourier space. 
Going to Fourier space and taking k small, we have 

L = i^^^^AiMX' + Ea^^s-fcs^) 2 



Dcfir 



mmg A ± = -^(A 1 ± A 
our Lagrangian becomes 



(51) 

- 1 K is diagonalized and 



Ana ^ 

9 

~ S 2 



■ h "(k,Al-k 3 A1) 2 



where m± = —3, —1 are the two eigenvalues of K. 
We can write this as 



L = J2(A^M q (A q ) 



(52) 



(53) 



where 



M q = 



( JL 

2g 3 a 
_ im q k2 

ira q k\ 
4-ira 
i kgk 3 
\ 2g 3 a 



im q k 2 
4-rra 
J 3 ak\ 
2 

im q ko 
4ira 
J 3 ak 1 k 3 
2 



im Q k 



qKl 



4ira 
im q ko 
4-Ka 
J 3 ak 3 
2 

J 3 ak 2 k 3 
2 



2g 3 a 



kpk 3 
2g 3 a 
J 3 ak 1 k 3 
2 

J 3 ak2k 3 
2 

_ J 3 a(kf + kl) 
2 



E 4& A h - A h + A Ui,i) 2 ( 4? ) Choosing the temporal gauge, A 



(54) 

0, we can reduce 



M q to the 3x3 submatrix 



where 



-2 -1 
-1 -2 



(48) 



M q = 



J 3 ah $ 


iniqko 


2 


4iva 


im q k() 


J 3 ak 3 


4ira 


2 


J 3 akik3 


J3ak 2 k 3 



^0 

2g 3 a 



J 3 aki_k 3 
2 

J 3 ak2k 3 
2 

_ J 3 a(kl+kl) 
2 



(55) 
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Setting the determinant to we find one gapless mode 
for each q = ± with dispersion 

= J 3 g 3 a\kl + kl) + A -^~ki (56) 

In principle, these gapless modes should be visible in in- 
elastic light scattering measurements. Thus, such mea- 
surements could be used to distinguish the 3D parton 
state from other candidate states - such as the decou- 
pled Laughlin state. In section [X] we describe another 
experimental signature of the 3D parton state, involving 
surface excitations. 

In addition to the gapless gauge excitations, the 3D 
state also contains gapped parton excitations. These ex- 
citations are fermions, and carry electric charge q = e/3 
as in the A^-layer case (|37[) . They are minimally cou- 
pled to the gapless gauge bosons and therefore have long 
range interactions. Note that these excitations are truly 
3D quasiparticles - they can propagate freely both within 
layers and between layers. This should be contrasted 
with the charge e/3 particles in the decoupled Laughlin 
state, which are confined to individual layers. In this 
sense, the parton state is "more 3D" than the decoupled 
Laughlin state. 

VIII. CROSSOVER FROM 2D TO 3D 

In the previous sections, we referred to the 3D limit as 
the limit N — > oo. However, the 3D limit can also be ac- 
cessed when there are a finite number of layers - provided 
that we probe the system at appropriate length and en- 
ergy scales. In this section, we discuss this crossover from 
2D to 3D physics. 

Consider an TV-layer system with N ^> 1. This system 
is described by three species of partons minimally coupled 
to the the gauge theory (|47|) . According to the analysis 
following (I47|) . the low energy modes of this gauge theory 
satisfy the dispersion relation (|56|) . Since N is finite, k 3 
is quantized in multiples of ir/Na. For each value of k 3l 
there is a corresponding 2D mode. 

The mode with the smallest gap corresponds to k 3 = 
ir/Na; the dispersion relation for this mode is 

- 2 -W(fc? + fcI) + ^ (57) 

We see that this mode has a gap A ~ J3/N 2 and a 
correlation length £ ~ N 2 a,yJ g 3 / J 3 . 

The gap A and correlation length £ are the impor- 
tant energy and length scales in the 2D/3D crossover. If 
one probes the system at energies less than A or lengths 
larger then £ (parallel to the layers), all the modes with 
&3 7^ will freeze out and the system will behave like 
a gapped 2D system. The physics is then described by 
the gapped FQH state in section IvTl On the other hand, 
if one probes the system at energies greater than A or 
lengths smaller then £, the system will behave like a 3D 
system. In this case, the physics is described by partons 
coupled to the gapless gauge theory (|47l) . 




FIG. 1: If r 3> £, the Berry phase associated with braiding a 
parton p around another parton p' is given by the 2D formula 
9 pp , = 2n/3N. If r <C £, it depends on the details of the path 
and can be calculated from the 3D gapless gauge theory (I47|l . 

One example of this crossover is the following thought 
experiment. Imagine adiabatically braiding one charge 
e/3 parton excitation p around another e/3 parton ex- 
citation p' - say of a different species - using a braiding 
path parallel to the layers (see Fig. [T|). First, consider 
the case where the separation r between the partons is 
kept larger than £. In this case, the modes ([5tl)) will be 
effectively frozen at this distance. While the presence of 
parton p' will change the gauge flux seen by parton p, 
the gauge flux will be localized to within a distance £ of 
p', and will be exponentially suppressed near the braid- 
ing path. Thus, the only interaction between the two 
partons will be a statistical interaction: the presence of 
parton p' will change the total gauge flux enclosed by the 
braiding path of parton p. The Berry phase associated 
with braiding one parton around the other is then given 
by the 2D mutual statistics formula: B pp i — 2ir/3N. 

Now, consider the case where the separation r between 
the partons is kept much smaller than £. In this case, the 
modes (|56[) will not be frozen out and the partons will 
experience long range interactions. The phase associated 
with braiding one parton around the other will depend 
on the details of the path, and can be calculated using 
the 3D gapless gauge theory (|4"T)) . 

As for the crossover between the two regimes, we ex- 
pect that the Berry phase in the 3D regime scales with 
the separation between the partons according to some 
power law. When the separation is of order £, we expect 
that the phase is of order 1 /N so that it agrees with the 
phase in the 2D regime. 

IX. RELATIONSHIP TO DECOUPLED v = 1/3 
LAYERS 

An interesting feature of the parton construction is 
that it can describe the decoupled v — 1/3 layered state 
within the same framework as the 3D parton state. One 
can tune from one state to the other by changing a single 
coupling constant in the parton gauge theory. 

To see this, let us go back to the original parton 
Hamiltonian ([5]). So far we have analyzed the physics 
of this Hamiltonian in the limit of weak gauge fluctua- 
tions: g 3 <C J3,t p3 ,t p ± and g± <C J±,t p3 ,t p ±. We found 
that in this regime, the low energy physics was described 
by the 3D parton state. 
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However, by increasing g$ one can also access a regime 
where (interlayer) gauge fluctuations are strong. More 
specifically, suppose that 53 3> J3,t p s,t p ±. In this case, 
the compactness of the U(l) gauge field becomes impor- 
tant: since the lattice electric field is integer valued and 
53 is large, the interlayer field E\ is essentially fixed at 
.Ef = 0. Nonzero values of E\ cost energy of order 33. 
As a result, the interlayer tunneling terms Ht z {z+\) an d 
interlayer flux terms cos(AiA* z 3 - A q xz i + A q , +1 , J in 
(|8|) are suppressed, and can be dropped from the Hamil- 
tonian. At low energies, the physics is then described by 
H = ~}2 z {H z a + H zt ) - the effective theory for decoupled 
v = 1/3 states. 

On an intuitive level, the basic physics is that large 
interlayer gauge fluctuations prohibit partons from tun- 
neling between the layers, leading to a decoupled layer 
state. One can also think about the transition between 
the 3D state and the decoupled layer state in terms of 
Higgs condensation. Consider, for example, the bilayer 
case. Recall that the interlayer gauge field terms ([23")) in 
the bilayer parton Hamiltonian are given by 

= £f(^V 3 ) 2 

xq 

- Y, J 3 COS ( A i A ll,3- A ll,i + A Ui) ( 58 ) 
xiq 

Let us view the operator e lj4 ' 13 as the creation operator 
of a boson at site x, while E xl 3 is the number operator 
which measures the number of bosons at site x. The first 
term is then a potential energy term which describes the 
energy associated with having a certain number of bosons 
on a given site, while the second term is a kinetic energy 
term which describes a boson hopping from one site to a 
neighboring site. The presence of the combination A\ i — 
A q 2 i in the argument of the cosine tells us that the boson 
is minimally coupled to the gauge field A\ i — A\ i . 

It is illuminating to think about the strong and weak 
gauge fluctuation regimes in this language. When 33 3> 
J3, the boson is massive (e.g. in a Mott insulating phase) 
and is therefore irrelevant at low energies. The interlayer 
gauge terms can then be dropped and the low energy ef- 
fective Hamiltonian (|19[) consists of two decoupled layers, 
H = 'Yli Z —x{H z t + Hza)- The ground state is thus two 
decoupled v = 1/3 states. 

On the other hand, when <C J3, the boson con- 
denses. Since the boson is minimally coupled to the 
gauge field A\ i — A\^ this boson condensation is a kind 
of Higgs condensation where the Higgs boson is coupled 
to a Chern-Simons gauge field. [H| When such a Higgs 
boson condenses, the result is another gapped FQH state 
- in this case, the bilayer parton state. 

Because one can tune from the decoupled state to the 
parton state by changing a single coupling constant 173, 
one can speculate that these states are in some sense 
neighboring or proximate phases. This is one of the rea- 
sons that we propose the parton state as a candidate for 
an intermediate tunneling phase. 



X. EDGE (AND SURFACE) STATES 

In this section, we discuss the edge states for the TV- 
layer system, both for finite N and in the 3D limit N — > 
00. 

First, consider the case of finite N. In this case, the 
edge theory can be read off from the AT-matrix and charge 
vector t described in section IVll using the standard for- 
malism. (27} The result is a chiral boson theory with 
3A — 2 modes, (j) 1 ,1 — 1, ...,3N — 2. The Lagrangian is 
of the form 

L = ^-{K u d t ^d x ^ - Vud^d^) (59) 

47T 

where V/j is a positive definite velocity matrix which 
describes the velocities of each of the modes and the 
density-density interactions between different modes. 
Quasiparticle excitations are parameterized by integer 
vectors I and are created by operators of the form 
exp(iZ/0 / ). The electric charge corresponding to a quasi- 
particle I is given by q = e ■ (t T A _1 Z). 

As many aspects of the edge theory depend on mi- 
croscopic details of the edge and can be affected by edge 
reconstruction, let us discuss two simple quantities which 
are universal. The first quantity - the electric Hall con- 
ductance - is given by 

(60) 

Of course, this is exactly what we expect since the parton 
state has v = 1/3 per layer. 

A more interesting quantity is the thermal Hall con- 
ductance. Recall that each chiral boson edge mode gives 

TT 2 k 2 

a contribution of ± „, B T to the thermal Hall conduc- 

3h 

tance, with the sign determined by the chirality of the 
mode. Thus, the thermal Hall conductance can be com- 
puted by counting the number of edge modes. Since the 
AT-matrix has 37V — 2 positive eigenvalues and no negative 
eigenvalues, there are 3A — 2 modes propagating in one 
direction and no modes propagating in the opposite di- 
rection. We conclude that, the thermal Hall conductance 
is 3 A - 2 (in units of K£f-T). 

Before concluding this section, let us briefly discuss 
the 3D limit, N — > 00. In this case, the boundary is 
two-dimensional, and the edge states are actually sur- 
face states. Let us focus on the most interesting kind 
of boundary: a boundary in the xz plane (with layers 
oriented in the xy plane). 

The analysis of the surface states is complicated by the 
fact that the bulk has gapless modes. Because of this, we 
will not analyze the surface in detail, but rather sketch 
the basic qualitative picture which is evident in mean 
field theory. In mean field theory, the surface modes are 
given by 3 species of non-interacting 2D fermions, which 
are chiral in the x direction but non-chiral in the z direc- 
tion. The modes form a "sheath" of chiral fermions with 
Fermi surface k x ~ t z cos(k z a). [l| When one goes be- 
yond mean field theory and includes gauge fluctuations, 
these fermions will become minimally coupled to the bulk 
"photon" mode. The gauge fluctuations will certainly 
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affect the surface theory; however, we know that they 
cannot gap out the surface modes entirely. Indeed, the 
gaplessness of the edge modes is protected by the nonzero 
electric Hall conductivity in the bulk (e 2 /3/i per layer). 
The surface sheath therefore resembles a chiral analog 
of the Halperin-Lee-Read state [28| which is protected 
against gap-forming instabilities by the topological char- 
acter of the bulk phase. 

These surface states may provide the simplest experi- 
mental signature of the 3D parton states. In particular, 
consider the z-axis surface longitudinal conductance a zz . 
In mean field theory, a zz behaves just like the conduc- 
tance of the layered integer quantum Hall system studied 
in [l[. Thus, a zz ~ const, as T — > 0. Including gauge 
fluctuations, we expect that this constant will be renor- 
malized, but a zz will remain finite at zero temperature. 
This should be contrasted with the behavior of a zz in the 
decoupled v = 1/3 Laughlin state. In that case a zz ~ T 3 
as T —> 0. 1] Thus, a measurement of a zz could, in prin- 
ciple, distinguish the 31? parton state from the decoupled 
Laughlin state. 

XI. CONCLUSION 

In this paper, we have constructed a candidate state 
for a multilayer FQH system with average filling v — 1/3 
per layer. We have proposed that the state may be re- 
alized in the "intermediate tunneling regime" where the 
interlayer tunneling strength is of the same order as the 
Coulomb energy e 2 Our construction is based on a 
slave particle approach, known as the "parton construc- 
tion." 

We have analyzed the state for both a finite number 
of layers N, and in the 3D limit, N — ► oo. In the case of 
a finite number of layers N, the state is a gapped FQH 
state and is described by a (3iV — 2) x (3iV — 2) K- matrix. 
Its quasiparticle excitations are anyonic and carry charge 
e/3. 

In the 3D limit, N — » oo, the state is more unusual. It 
supports two types of excitations: gapped e/3 fermionic 
quasiparticle excitations and gapless neutral collective 
modes. The quasiparticle excitations are truly 3D quasi- 
particles and can propagate freely both within and be- 
tween layers. (This is in contrast to the charge e/3 exci- 
tations in the multilayer Laughlin state, which are con- 
fined to individual layers) . The gapless neutral collective 
modes are emergent "photon" modes originating from 



the slave particle gauge theory. Unlike Maxwell pho- 
tons, they come in only one polarization, and have an 
anisotropic dispersion w 2 . ~ + k z . The e/3 fermionic 
quasiparticles are minimally coupled to these "photon" 
modes so that they have long range interactions. 

The edge physics of the finite layer and 3D systems 
is also interesting. When N is finite, the edge theory is 
described by a conformal field theory with 3N — 2 chiral 
boson modes. In the 3D limit, the edge modes are more 
complex. In mean field theory, the edge (or more ac- 
curately, surface) modes are described by three different 
species of non-interacting 2D fermions which propagate 
chirally in the rr-direction (e.g. the direction parallel to 
the layers) and non-chirally in the z direction (e.g. the di- 
rection perpendicular to the layers). Going beyond mean 
field theory, we expect that these fermions are minimally 
coupled to the gapless bulk "photon" modes. However, 
we have not analyzed the surface physics in detail. This 
is an interesting direction for future research. 

Another direction for future research would be to con- 
struct other types of layered FQH states. For example, it 
would be interesting to build multilayer states with aver- 
age filling v = 1/2 per layer - in particular, multilayer 
states which are related to the Moore- Read v = 1/2 
state [33| or the composite Fermi liquid v = 1/2 state 
(28| (instead of the Laughlin v — 1/3 state which we 
have investigated here). One possible approach for this 
problem would be to employ a parton construction where 
one writes an electron as c = did^f, where d\^2 are 
fermionic partons carrying charge e /2 and / is a neutral 
fermionic parton. One could then consider mean field 
states where, in each layer, di, d2 are in integer quantum 
Hall states and / is in a p + ip superconducting state or 
a Fermi liquid state. In this way, one may be able to 
construct multilayer and 3D relatives of the Moore-Read 
or composite Fermi liquid states. 

In general, there are clearly many possibilities for 3D 
multilayer FQH states, most of which have not been ex- 
plored. We hope that the parton construction provides 
a useful tool for constructing and analyzing these new 
states of matter. 
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